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Abstract:  A  classical  force  field  has  been  developed  to  represent  die 
conformational  characteristics  of  the  polycarbonate  of  4,4'- 
isopropylidenediphenol  ("bisphenol-A  polycarbonate".  PC)  based  upon  recent 
experimental  and  quantum  mechanical  data  reported  for  dipbenylcarbonate 
and  diphenylpropane.  This  force  field  is  an  improvement  upon  previously 
published  molecular  mechanics  force  fields  because  it  allows  for  rotation 
about  all  the  single  bonds  in  the  PC  repeat  unit.  A  rotational  isomeric  state 
model  of  PC  has  been  obtained  using  the  force  field  results;  the  computed 
unperturbed  chain  dimensions  of  <r^>o/M  =  1.03  agree  well  with  reported 
experimental  values. 
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I.  INTRODUCTION  AND  PREVIOUS  WORK 


In  or3eT  to  simulate  dense  structures  of  poly (oxycarbonyloxy- 1,4- 
phenylenc-isopropylidene-l,4-phenylene)  ("bisphcnol-A  polycarbonate",  PC), 
an  essential  requirement  is  detailed  knowledge  of  both  the  short  range  and 
long  range  interactions  in  the  system.  In  the  following  we  develop  the  atom- 
atom  force  field  necessary  and  explore  the  effects  of  the  short  range 
(intramolecular)  interactions  by  investigating  the  conformational 
characteristics  of  the  fragment  molecules  diphenyl  carbonate  (DPC)  and  2,2- 
diphenylpropane  (DPP).  These  molecules  were  chosen  to  represent  the 
conformational  characteristics  of  the  PC  repeat  unit.  In  a  companion  paper  to 
this  one  we  report  on  a  simulation  of  the  dense  structure  of  PC  performed  in 
collaboration  with  Gentile  and  Ludovice  which  in  the  remainder  of  the 
present  paper  will  be  referred  to  as  II.* 

The  conformational  characteristics  of  DPC  and  DPP  have  been  the 
subject  of  several  detailed  studies  involving  molecular  mechanics^'^  and 
quantum  mechanics  methods.***'*^  Previous  work  on  molecular  mechanics 
could  not  benefit  from  the  recent,  detailed  experimental  results  and  quantum 
calculations,  and  therefore  does  not  accurately  describe  the  conformations  of 
the  fragments  of  PC.  One  of  the  earliest  works  published  on  the  conformational 
analysis  of  DPC  and  DPP  using  molecular  mechanics  is  that  by  Williams  and 
Flory.^  Their  approach  is  based  on  structural  data  from  tables  of  interatomic 
distances*^  and  a  6-12  power-law  potential  energy  function  for  the  nonbonded 
interaction!  between  pairs  of  atoms.  Despite  the  simplicity  of  the  approach, 
the  resulting  conformational  analysis  was  quite  accurate;  however,  rotational 
barriers  and  electrostatic  interactions  were  neglected  in  these  calculations. 
The  later  work  by  Tonelli,^  Erman  et  al..^>^  and  Sundararajan^  is  based  on  the 
<‘x^)enmental  aata,  with  additional  information  on  DPC  from  crystal 


2 


structure  analysis. 20  The  nonbonded  interactions  used  in  these  calculations 
have  employed  a  Lennard- Jones  6-12  and  Hill's  functions.  Although  these 
studies  go  into  more  detail  than  the  earliest  work,  such  as  including  some 
single  bond  torsional  energy  functions,  no  electrostatic  interactions  are 
included,  and  for  DPC  the  carbonate  group  is  considered  rigidly  fixed  in  the  all 
trans  conformation.  Recent  experimental  and  theoretical  data  clearly  indicates 
that  this  assumption  is  invalid.  Tckely  and  Turska^  have  attempted  to  rectify 
this  problem,  but  at  the  time  did  not  have  available  sufficient  data  to  obtain 
values  for  all  the  relevant  force  field  parameters. 

Recently  Perez  and  Scaringe  presented  results  utilizing  single  crystal  X- 
ray  diffraction  data  on  the  dicarbonate  4,4'-isopropylidenediphenylbis- 
(phenylcarbonate)  (DBPC).*  Henrichs,  Luss  and  Scaringe  have  further  studied 
the  crystalline  structure  of  DBPC,  identifying  a  second  crystal  structure.^  The 
results  of  the  two  studies  are  very  informative. 

Considerable  research  has  been  reported  applying  quantum  mechanics 
to  the  conformational  and  structural  analysis  of  DPC  and  DPP.  Among  the 
several  semi-empirical  methods  used  on  DPC  and  DPP  are  CNDO  (Complete 
Neglect  of  Differential  Overlap).  MNDO  (Modified  Neglect  of  Differential 
Overlap),  and  AMI  (Austin  Model  1).*®’*'^  Bicerano  and  Clark^^*^^  have 
published  two  very  detailed  studies  which  include  not  only  results  obtained 
with  MNDO  and  AMI,  but  also  PRDDO  (Partial  Retention  of  Diatomic  Differential 
Overlap)  calculations,  an  approximate  ab  initio  method,  on  DPC,  DPP,  and  on 
the  dicaiixMute  of  4,4'-isopropylidenediphenol  and  two  methanols.  Of  the  ah 
initio  methods,  STO-3G  level2l>22  Hartree-Fock  calculations  hi.e  been  the  most 
widely  employed. Laskowski  et  al.^^  and  Jaffe^^  have  also  reported  the 
results  of  calculations  on  DPC,  DPP,  and  smaller  model  compounds  at  the  STO-3G, 
4-31G.21.23  and  6-31G*  levels.21.23 
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On  the  whole,  there  exists  consistent  experimental  data  and  reliable 
theoretical  calculations  on  DPC  and  DPP.  On  the  basis  of  these  studies,  a  new 
molecular  mechanics  method  can  now  be  developed. 

II.  ENERGY  REPRESENTATION 

The  molecular  structure,  i.c.  the  bond  lengths  and  bond  angles,  was 
chosen  based  upon  the  experimental  and  theoretical  data  available  in  the 
literature.  A  Lennard-Jones  potential  energy  function  was  selected  to 
represent  the  van  der  Waals  interactions^^  and  a  Coulombic  potential  with  a 
distance  dependent  dielectric  constant  was  used  to  represent  the  electrostatic 
interactions. 25  Structure  specific  intrinsic  torsional  potential  energy 
functions  have  been  included  to  represent  the  bonded  interactions  (see 
below).  Since  the  purpose  of  this  force  field  was  primarily  for  its  use  in  the 
generation  of  dense,  glassy  microstructures  of  PC  (described  in  II),  the 
assumptions  of  fixed  bond  angles  and  bond  lengths  and  the  representation  of 
the  methyl  groups  as  spherical  pseudo^atoms  were  made  for  computational 
expedience.  The  choice  of  the  intrinsic  torsional  energy  functions  is  strongly 
dependent  on  these  assumptions. 

Figure  1  shows  the  fragment  from  a  PC  structure  in  the  conformation 
where  all  torsion  angles  are  zero.  The  geometric  parameters  used  are  described 
in  Figure  I  and  in  Table  I.  They  are  based  upon  the  X-ray  diffraction 
experimenta  on  DBPC  by  Perez  and  Scaringc*  and  by  Henrichs  et  al..^  the 
diffraction  experiments  by  Yoon  and  Flory  on  DPC.^O  and  the  quantum 
mechanical  calculations  of  Bicerano  and  Clark.  ^5  Because  the  model  developed 
here  does  not  allow  for  any  hond  or  changes  in  the 

length,  some  adjustments  needed  to  be  made  if  realistic  torsional  barriers  were 
to  be  obtained;  The  bond  angles  of  the  carbonate  group  were  chosen  to  be 


slightly  larger  than  the  experimental  values.  The  numbering  scheme  for  the 
bonds  and  torsion  angles  described  in  Figure  1  will  be  used  throughout  this 
paper. 

The  nonbonded  forces  include  van  der  Waals  and  electrostatic 
interactions.  The  parameters  for  the  nonbonded  potential  energy  functions 
representing  these  interactions  are  listed  in  Table  II. 

A  Lennard-Jones  function  is  used  to  represent  the  van  der  Waals 
interactions: 

where  is  the  pairwise  van  der  Waals  interaction  between  atoms  i  and  j,  and 

O  is  the  fundamental  length  measure  of  the  potential  where  the  interaction 

vanishes  and  is  obtained  using  the  Slatcr-Kirkwood  equation.^* 

The  electrostatic  nonbonded  interactions  arc  represented  by  the 
following  Coulomb  potential: 


where  U^jj  is  the  electrostatic  pairwise  interaction  potential,  E©  is  the 
permittivity  in  vacuum,  D(r)  is  the  effective  dielectric  constant,  and  the  q's 
are  the  partial  charges  on  atoms  i  and  j.  At  atomic  level  dimensions,  the 
dielectric  constant  is  strongly  dependent  on  the  interatomic  separation,  and 
therefore  a  distance-dependent  dielectric  constant  D(r)  is  used,  based  on  an 
approximation  first  introduced  by  Block  and  Walkcr:^^ 


CKr)  =  1 


for  r  <  a 


(^a) 

(3b) 


DTf)  =  £b  for  r  >  a 

K  =  aln(eB)  (3c) 

where  £8  is  the  dielectric  constant  of  the  bulk  polymer,  which  for  these 
calculations  was  taken  to  be  3.07.^®  The  value  a  of  r  at  which  D  switches  from 
unity  to  an  exponentially  decaying  function  was  chosen  to  be  the  same  for  all 
atom  pairs,  i.  e.,  a  =  3.0S  A.  This  value  is  the  sum  of  the  smallest  and  the  largest 
van  der  Waals  radii  for  the  charged  atoms.  The  choice  of  a  made  here  differs 
from  that  made  in  the  calculations  of  Ludovice  and  Suter,^^  where  a  was 
allowed  to  vary  according  to  the  sum  of  the  van  der  Waals  radii  of  atoms  i  and  j. 
The  coulombic  potential  including  the  Block-Walker  approximation  for  two 
hydrogens  is  shown  in  Figure  2.  The  partial  point  charges,  qi,  shown  in  Table 
II.  are  based  upon  MNDO  and  AMI  calculations,^^  and  were  chosen  so  that  the 
PC  repeat  unit  is  electrically  neutral.  A  rough  test  for  the  partial  charges  can 
be  obtained  by  comparing  computed  and  experimental  values  for  the 
considerable  dipole  moment  of  DPC.  The  estimated  value  from  the  force  field 
(keeping  DPC  electrically  neutral)  is  0.97  Debye,  and  shows  good  agreement 
with  the  experimental  value  of  0.87  +  0.05  Debye.'*  Similarly,  the  charges  in  the 
phenylene  rings  can  be  used  to  obtain  an  estimate  of  the  quadrupole  moment 
of  benzene;  the  calculated  value  is  -29.0  x  10*'**^  Coulombs  m^  which  compares 
favorably  with  the  experimental  value  of  -(29.0  ±  1.7)  x  10"'*®  Coulombs  m2.3  2 
In  addition  to  the  nonbonded  interactions,  intrinsic  torsional  potentials 
must  be  included  to  complete  the  description  of  the  interactions  along  the 
molecule.  Such  bonded  interactions  are  typically  represented  by 
trigonometric  functions  and  are  chosen  so  that  the  bond  rotation  barriers  and 
the  conformational  minima  determined  with  the  force  field  closely  resemble 
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experimental  or  theoretical  values.  In  the  cases  of  DPC  and  DPP,  the 
trigonometric  "functions  chosen  were  of  somewhat  unusual  form. 

For  DPC,  the  intrinsic  torsional  energy  function  about  the  Car  .  o  bond 
(a  carbonate  flanking  bond)  is  a  simple  squared  sinusoidal  function,  often 
used  to  represent  the  effect  of  ic -orbital  overlap  on  conformation,  and  is  shown 
in  Figure  3a: 


U<p,(  <p)=  B<p,sin2(p 


(4) 


d<p 


=  B<,>,  sin  2  <p 


(5) 


where  U(p^  is  the  intrinsic  torsional  energy.  is  the  barrier  height,  and  <p  is 

the  rotational  angle  measured  from  the  datum  plane  of  the  molecule  as 
depicted  in  Figure  1. 


For  rotation  in  the  carbonate  group  about  the  0  -  CC*0  bond,  it  was 
necessary  to  devise  a  more  complex  function,  implicit  in  the  dummy  variable  x: 

9  =  ^  (  T  -  X  sin  x) 

U<p4=^(l  -  cost) 
sin  x 

( 1  -  X  cos  x) 

where  U(p^  is  the  intrinsic  torsional  energy  ,  B(p^  is  again  the  barrier  height, 
and  9  is  the  rotational  angle.  The  parameter  X  controls  the  width  of  the  energy 
well  (if  X  =  0  the  function  is  a  simple  sinusoidal,  and  if  X  =  1.  the  function  is  a 


(6a) 

(6b) 

(7) 
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cycloid).  Here,  X  was  chosen  to  be  0.3.  The  intrinsic  torsional  energy  function 
is  shown  in  Frgnre  3b.  The  need  for  such  a  torsional  energy  function  arose 
from  the  desire  to  keep  the  minimum  energy  conformaiiouo  planar  while  at 
the  same  time  keeping  the  total  rotational  energy  barrier  about  the  bond  low 
(see  below). 

For  the  rotation  about  the  Car  .  cal  bonds  (the  bonds  joined  to  the 
isopropylidene  group)  it  was  again  necessary  to  devise  a  particular  functional 
form  for  the  energy  minima  to  be  correctly  located: 


U, 


=?^ri 
<P1  2 


cos  {  4  { <p  -  V  sin  2(p  ) )  ] 


(8) 


=  -2  B<p,  sin  {  4  ( <p  -  V  sin  2((>  ) )  ( 1  -  2  v  cos  2(p ) 

where  Uepj  is  the  intrinsic  torsional  energy,  B(p  j  is  the  barrier  height,  and  V  is 

a  parameter  which  changes  the  location  of  the  energy  minima.  This  function, 
where  v  was  set  to  0.3,  is  shown  in  Figure  3c. 

III.  CONFORMATIONAL  ANALYSIS  OF  THE  MODEL  COMPOUNDS 
3.1  DPC 

DPC  has  two  sets  of  conformations  of  low  energy,  the  more  favored  one 
occurring  when  the  O  -  CC=0  bonds  of  the  carbonate  group  are  in  their 
trans  ,trans  conformation.  At  the  secondary  set  of  minima  the  carbonate  group 
is  in  its  cis, trans  or  its  trans, cis  conformation  and  has  a  slightly  higher 
energy.  The  carbonate  group  itself  is  planar  due  to  electron  delocalization.  The 
phenyl  rings  are  rotated  out  of  the  plane  of  the  carbonate  group  because  of 
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opposing  forces.  Electron  delocalization  between  the  carbonate  group  and  the 
phenyl  rings'  favor  an  all  planar  conformation,  whereas  the  steric  hindrance 
between  the  ortho-hydrogen  on  the  phenyl  ring  and  the  carbonate  oxygen 
favor  a  perpendicular  arrangement  of  the  phenyl  rings  and  the  carbonate 
group.  The  rotation  of  the  -  O  bonds  arc  not  coupled;  all  the  combinations  of 
cither  positive  rotations  or  negative  rotations  have  approximately  the  same 
energy. 

In  Table  III  the  force  field  results  developed  here  are  compared  with  the 
published  experimental  and  theoretical  values.  The  location  of  the  minima 
refer  to  the  dihedral  angles  of  the  phenyl  rings,  where  all  the  angles  are  zero 
in  the  all  planar  conformation  shown  in  Figure  1.  The  values  of  AE  =  Eds  - 

Etrans  kcal/mol  refer  to  the  difference  in  the  energy  between  the 

trans,trans  conformation  and  the  cis.trans  or  the  trans.cis  conformation.  The 
achieved  agreement  between  our  adjusted  force  field  results  and  the 
corresponding  theoretical  and  the  experimental  values  from  the  literature  are 
satisfactory. 

Furthermore,  the  energy  barriers  to  rotation  were  considered.  In  Table 
IV,  the  values  calculated  from  the  developed  force  field  are  compared  with 
available  theoretical  values.  The  bond  energy  barrier  for  the  -  O  bond  is 
plotted  in  Figure  4  and  shows  two  energy  barriers  to  rotation:  the  smaller 
barrier  at  (90”, 2/0”)  is  usually  attributed  to  the  breaking  up  of  the  electron 
delocalization  between  the  phenyl  ring  and  the  carbonate  group,'^  while  the 
larger  barrier  at  (0°,180”)  is  due  to  the  steric  hindrance  between  the  ortho- 
hydrogen  and  the  carbonate  oxygen.  The  contour  map  shown  in  Figure  S 
displays  the  energy  as  a  function  of  the  two  -  C  bonds  (angles  94  and  95) 

in  the  carbonate  group.  The  map  clearly  illustrates  that  the  carbonate  group 
being  much  'softer'  than  previously  assumed,  must  be  represented  with  a 


9 


higher  degree  cf  flexibility  than  was  done  in  earlier  force  fields.  The  force 
field  results  agree  with  the  theoretical  values,  except  that  the  larger  barrier 
of  the  Ca*"  -  O  bond  is  slightly  higher  than  the  theoretical  results. 

3.2  DPP 

The  DPP  molecule  has  one  symmetrically  unique  conformational 
energy  minimum,  where  the  phenyl  rings  rotate  out  of  the  plane  in  a 
propeller-like  conformation.  This  rotation  is  due  to  the  steric  hindrance 
between  the  ortho-hydrogens  on  different  ohenyl  rings.  In  Table  V  the  force 
field  results  devised  here  are  compared  with  available  experimental  and 
theoretical  values.  The  minimum  refers  to  the  dihedral  angle  of  the  phenyl 
rings,  where  the  all  planar  conformation  is  defined  as  (0°,  O^i.  The  energy 
barrier  refers  to  the  bond  total  rotational  energy  barrier  to  phenyl  ring 
rotation  about  the  bond.  The  optimal  path  for  rotation  is  when  the 

phenyl  rings  move  simultaneously  through  a  saddle  point  at  (0®.90°)  or 
(90® .0°).  In  Figure  6,  the  force  field  results  (a)  are  compared  with  the  STO-3G 
level  results  of  Laskowski  et  al.^^  The  energetically  preferred  path  of  rotation 
is  indicated  by  the  dotted  line. 

IV.  ROTATIONAL  ISOMERIC  STATE  MODEL 

From  the  results  obtained  from  the  fine  tuned  force  field  described 
above,  a  Rotational  Isomeric  State  (RIS)  model  was  deduced. The  statistival 
weight  matrices  thus  obtained  are  shown  in  Table  VI,  where  the  bo.^d 
numbering  refers  to  the  scheme  used  in  Figure  1.  The  parameter  Y  gives  the 

fractional  concentration  of  the  cis.trans  or  trans.cis  states  as  a  function  of  AE, 
the  difference  in  energy  between  the  cis.trans  or  the  trans.cis  state  and  the 
energy  of  the  trans.trans  state,  AE  =  Ecij  -  Efrans  •  *5* 
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Y  =  e-p  '  -  AE/  RT) 


(10) 


(The  value  for  Y  considered  here  is  slightly  different  from  the  one  reported  by 
us  earlier^^)  The  cis.cis  conformation  brings  the  two  phcnylene  rings  into 
very  close  proximity,  and  the  concomitant  excess  energy  makes  this 
conformation  inaccessible. 

Previous  rcr  arch  that  utilized  the  RJS  method  for  PC  includes  work  by 
Williams  and  Flory^  and  by  Laskowski  et  al.*^  Both  groups  of  researchers 
employed  the  same  RIS  scheme,  a  model  with  significantly  fewer  states  and 
fewer  expl.*-itly  considered  bonds  than  the  one  presented  here.  Our  RIS 
scheme  differs  from  that  of  Williams  and  Flory  in  tha;  we  describe  the 
conformations  of  the  repeat  unit  with  6  statistical  weight  matrices,  whereas 
Williams  and  Flory  use  only  4.  Also,  we  include  a  4  state  model  which  explicitly 
represents  the  symmetrically  equivalent  positions  of  the  phenylenc  ring 
conformations  about  the  carbonate  group  and  about  the  isopropylidene  group. 

The  conformations  of  the  phenylene  rings  about  the  carbonate  group 
are  mutually  independent  and  have  a  4-fold  symmetry,  resulting  in  4 
equivalent  states  (see  U3,  in  table  VI).  All  first  order  interactions  arc 
identical  and  no  second  order  parameter  is  required. 

Two  matrices  are  used  to  represent  the  conformations  of  the  phenylene 
rings  about  the  isopropylidene  group.  The  first  matrix,  U  ] ,  encompasses  4 
equivalent  states  representing  the  rotational  minima  found  ii  DPP  which  are 
independent  of  the  previous  bond.  The  second  matrix,  U^,  comprises  4  states 
which  are  dependent  on  the  state  chosen  for  the  previous  bond  so  that  only 
(he  "p.opeiler-like”  conformations  of  the  rings  is  allowed  (i.  e.  ring  rotations 
of  (-*45°. +45'')  are  allowed,  whereas  ring  rotations  of  (+45°,+ 135°)  are  not). 
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Although  our  scheme  and  that  of  Williams  and  Flory  differ  formally, 
using  equivalent  conformational  parameters  the  two  schemes  yield  identical 
results.  The  RIS  scheme  presented  here,  even  though  it  is  of  lesser  elegance, 
has  the  major  advantage  that  it  can  be  easily  modified  to  allow  for  simple 
substitution,  extending  its  use  to  other  polycarbonates  and  to  other 
structurally  similar  polymers.  Our  RIS  scheme  representation  can  also  be 
used  without  modification  for  the  construction  of  initial  chain  conformations 
used  for  the  generation  of  dense,  amorphous  microstructures  as  we 
demonstrate  in  II. 

Values  calculated  for  the  limiting,  unperturbed  mean-square  end-to- 
end  distance  for  PC  by  standard  generator  matrix  techniques'^  are  compared 
with  experimental  values  in  Table  VII.  The  value  of  1.03  deduced  from  the 
above  described  force  field  compares  very  favorably  with  experimental  data. 

V.  CONCLUSIONS 

A  classical  force  field  has  been  developed  for  the  PC  repeat  unit  based 
upon  the  experimental  and  theoretical  data  that  has  recently  become 
available.  The  inclusion  of  electrostatic  charges  and  structure  specific 
torsional  potentials  creates  a  force  field  more  detailed  than  those  previously 
published.  The  principle  difference  between  this  force  field  and  earlier  ones 
is  that  torsion  around  all  single  bonds  is  considered.  Earlier  molecular 
mechanics  models  had  fixed  the  carbonate  group  in  its  all  planar 
conformatioa  or  did  not  allow  for  full  rotation  about  the  O  -  CC=0  bonds.^'*^’* 

The  results  show  that  there  is  good  agreement  between  the 
experimental  and  the  quantum  mechanical  data.  There  is  also  good  agreement 
between  results  obtained  with  this  force  field  and  quantum  mechanical  data 
with  respect  to  the  location  of  the  energy  minima  and  their  relative  energies, 
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and  also  with  respect  to  the  energy  barriers  to  rotation.  The  potential 
weakness  of  the  force  field  is  the  use  of  fixed  bond  angles,  fixed  bond  lengths, 
and  the  representation  of  the  methyl  groups  as  spherical  pseudo-atoms.  These 
adjustments  are  computational  expedients,  implemented  with  regards  to  the 
intended  use  of  this  force  field  in  dense  packing  simulation  (see  II).  and  have 
been  taken  into  consideration  when  the  geometrical  parameters  were  chosen. 
The  intrinsic  torsion  potentials  were  selected  to  compensate  as  far  as  possible 
for  these  shortcomings. 

A  new.  and  very  flexible,  RIS  model  was  thus  deduced  from  the  results  of 
the  force  field  computations.  Computed  values  for  the  characteristic  ratio  of 
the  unperturbed  chain  dimensions  compare  very  well  with  experimental 
values  in  the  literature. 
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Data 


Table  I.  Geometrical  Data  of  the  PC  Repeat  Unit. 


of  Relevance  for  the  Rotational  Isomeric  State  Sctiemf 


bond  length 


ll 

1.54 

12 

1.54 

1* 

2.76 

13 

1.41 

14 

1.33 

15 

1.33 

l6 

1.41 

Bond  Angles  (deg) 


71-02 

109.5 

71-04 

124.0 

7t  -  05 

109.0 

7C  -  06 

124.0 

xAi 


Details  of  the  Molecular  GeomcllY 


bond  length  (A) 
QOT  .  car  1.38 

-  Methyl  1-53 

CC=0.oC=0  1.21 

Car.H  110 


Bond  Angles  (deg) 

Z  120.0 

^0-CC=0.oC=0  125.5 

Z  C“  -  cal  -  Methyl  109  0 
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Table  II.  Potential  Energy  Functions  Parameter  Values. 


!■ 

atom  or  group  * 

ri°  .b  (A) 

BIIEH 

IH 

1 

Cal 

5.0 

0.00 

Car,  H  substituent 

1.85 

1.23 

5.0 

-0.14 

Car,  c  substituent 

1.85 

1.23 

5.0 

0.00 

Car,  0  substituent 

1.85 

1.23 

5.0 

0.08 

CP=0 

1.70 

1.23 

5.0 

0.50 

H 

1.20 

0.42 

0.9 

0.14 

-0- 

1.50 

0.70 

7.0 

-0.18 

8 

Methyl 

1.90 

1.77 

7.0 

0.00 

9 

ooo 

1.50 

0.84 

7.0 

-0.30 

^  Cal  identifies  an  aliphatic  carbon,  an  aromatic  carbon;  the  superscript 
'C=0'  indicates  atoms  that  arc  part  of  the  carbonyl  group, 
b  rj®  =  Van  dcr  Waals  radius.26 
c  ai  =  Polarizability. 

Ne  =  Effective  number  of  electrons.^* 

^  qi  =  Partial  electrostatic  point  charge.^ ^ 
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Table  III. 

Comparison  of 

Energy  Minima  in  DPC. 

Method 

Principal  M 

iflilQA  (<P3.<P6>  deg)  A  E  .(kcal/moh 

X-ray, 

48,48 

- 

AMI,  13 

44,44 

- 

PRDDO,  13 

44,44 

1.1 

STO-3G.  l"^ 

57.57 

1.68 

6-31G*.  b 

- 

2.74,  (2.21) 

This  Force  Field 

47,47 

1.7 

*  AE  =  Eq\s  -  Etrans 

b  The  6-3 IG*  calculations  arc  from  reference  18  and  were  done  for 
methylphcnyl  carbonate  using  two  different  methods. 
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Method 


Table  IV.  Total  Rotational  Energy  Barriers  in  DPC. 

Barrier  .  Q  (kcal/mol^  ^  Barrier  Q  -  C<^=Q  (kcal/mol)  ^ 


PRDDO. 

0.6 

4.0 

STO-3G.  17 

1.23 

4.5 

This  Force  Field 

3.3. 

(0.86) 

4.0 

^  Barrier  for  the  rotation 

of  the 

car  -  0 

bond. 

all  other  bonds  being 

continuously  adjusted  to 

give 

minimum 

energy 

conformations. 

I’  Barrier  for  the  rotation 

of  the 

0  -  C9=o 

bond. 

all  other  bonds  being 

continuously  adjusted  to 

give 

minimum 

energy 

conformations. 
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Table  V.  Energy  Minima  and  Total  Rotational  Energy  Barrier  in  DPP. 


Method 

Minima.  «pi.92.<ieg) 

Barrier  Car  .  cal  (kcal/moll 

X-ray, 

55.55 

- 

AMI,  13 

53,53 

- 

PRDDO,  13 

48.48 

- 

STO-3G.  17 

51,51 

1.9 

This  Force  Field 

45,45 

2.7 

^  Barrier  for  the  rotation  of  the  -  C®*  bond,  all  other  bonds  being 
continuously  adjusted  to  give  minimum  energy  conformations. 
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Table  VI.  Rotational  Isomeric  State  Parameters  for  PC 


Bond  Number  i 

ai 

0 

{  (pj_l  .  decrees 

Hi  a 

1 

45.135.225.315 

n  1  I  ll 

1  1  1  1  11 

1  1  1  I  1  1 

Ll  1  1  iJ 

2 

70.5 

45.135.225.315 

r  1 0 1  oi 
10  1  0  1 1 

1  1  0  1  ol 
Lo  1  0  ij 

3 

0 

45.135.225.315 

n  1 1  ll 

1 1 1 1  ll 

1 1 1 1  ll 
Ll  1  1  iJ 

4 

56.0 

0.  180 

ri  y1 

1 1  yI 

1 1  yI 

Ll  yJ 

5 

71.0 

0.  180 

1 — 

O 

6 

56.0 

45.135.225.315 

ri  1  1  ll 
Ll  1  1  iJ 

2  Y  =  exp  (-900A')  with  T  in  K,  at  room  temperature.  Y  “  0.05. 
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Table  VII,  The  Limiting,  Unperturbed  Mean-square  End-to-end  Distance  of 

PC. 


Method  lim  <  r  ^  M  (A  mol  g 

M  — >  oo 

This  RIS  Model  1.0? 


Light  scattering 


0.75  a 
0.98  b 


[  n  ]  0.87  a 

SANS  1.25  c 

1.28  d 


^  Reference  37. 
^  Reference  38. 

Reference  39. 
^  Reference  40. 
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ss 


(loin/t^oji)  XSjsua 


Figure  2.  The  potential  energy  of  electrostatic  interactions  for  two  hydrogen  atoms,  and 
the  Block  and  Walker  interpolation  approximation. 


Figure  3.  Dependence  of  intrinsic  torsional  potential  energy  functions  on  specific  torsion 
angles,  a),  dependence  of  U(p3  on  93  corresponding  to  the  C*'’  -  0  bond  in  DPC, 

b).  dependence  of  U(p4  on  94  referring  to  the  O-CC-0  bond  on  DPC.  and  c).  dependence 
of  U(p^  on  (pi  describing  the  situation  for  the  Caf  -  C*!  bond  in  DPP. 


3.3 


(lom/po^i)  ASaaua 


Figure  5.  Calculated  total  rotational  energy  for  DPC  as  a  function  of  the  torsion  angles  94 
and  (p5  labeled  in  Figure  1 .  The  values  of  all  other  torsion  angles  are  continuously 
adjusted  to  give  minimum  energy  conformations.  Contours  are  drawn  with  respect  to  the 
absolute  minima  at  94 . 95 . 0*  {trans.trans  conformation)  with  intervals  of  1 
kcal/mol.  X  denotes  minima,  and  ^  denotes  cols  between  troughs. 


